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Abstract 

The production of beauty quarks with a D*^ and a muon in the final state has 
been measured with the ZEUS detector at HERA using an integrated luminosity 
of 114 pb~^. Low transverse-momentum thresholds for the muon and D* meson 
allow a measurement of beauty production closer to the production threshold 
than previous measurements. The beauty signal was extracted using the charge 
correlations and angular distributions of the muon with respect to the D* meson. 
Cross sections for photoproduction and deep inelastic scattering are somewhat 
higher than, but compatible with, next-to-leading-order QCD predictions, and 
compatible with other measurements. 
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1 Introduction 



The production of beauty quarks in ep collisions at HERA is a stringent test for per- 
turbative Quantum Chromodynamics (QCD) since the large 6-quark mass {nib ~ 5 GeV) 
provides a hard scale that should ensure reliable predictions. For 6-quark transverse mo- 
menta comparable to the 6-quark mass, next-to-leading-order (NLO) QCD calculations 
in which the b quarks are generated dynamically are expected to provide accurate predic- 
tions [1,2]. 

The beauty-production cross section has been measured in pp collisions at the SppS [3] 
and Tevatron colliders [4], in 77 interactions at LEP [5], in fixed-target nN [6] and pN [7] 
experiments, and in ep collisions at HERA [8-12]. While most results, including recent 
results from the Tevatron, are in agreement with QCD predictions, some, in particular 
those from LEP, show large discrepancies. 

This paper reports a measurement of beauty production via the reaction ep ebbX — >■ 
eD*fiX' using the ZEUS detector at HERA. This reaction offers the advantage of providing 
a data sample enriched in b quarks and with strongly suppressed backgrounds from other 
processes, which allows low-p^ threshold cuts to be applied. This analysis therefore yields 
a measurement of beauty production closer to the production threshold than previous 
HERA measurements based on leptons and/or jets with high transverse momentum [8-11]. 
A similar measurement has been performed by the HI collaboration [12]. 

Of particular interest are events in which the muon and D* originate from the same parent 
B meson (Fig. [T^), e.g. —>■ D*~ ^'^u^. These yield unlike-sign D*-muon pairs produced 
in the same detector hemisphere. Due to the partial reconstruction (e.g. missing neutrino) 
the invariant mass is constrained to lie below the 5-meson mass. Another important 
contribution arises from charm-pair production, where one charm quark fragments into a 
D* and the other decays into a muon (Fig. Wp)- This again yields unlike-sign D*-muon 
pairs, but with the D* and the muon produced predominantly in opposite hemispheres. 
In addition, beauty-pair production in which the D* and muon originate from different 
beauty quarks can yield both like- and unlike-sign D*-muon combinations, depending on 
whether the muon is from the decay of the primary beauty quark (Fig. [It), or from a 
secondary charm quark (Fig. [T]i), and whether B^-B^ mixing has occured. 

Background contributions to both like- and unlike-sign combinations include events with 
either fake D* mesons, originating from combinatorial background, or fake muons. In this 
analysis, the signal is extracted from the unlike-sign sample, while the like-sign sample is 
used as a cross check. 

Cross sections are extracted separately for the photoproduction (7/?, photon virtuality 
< 1 GeV^), and deep inelastic scattering (DIS, > 1 GeV^) regimes, as well as 
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for the entire range in Q^, which includes the kinematic region in which these two event 
classes cannot easily be distinguished. 

2 Experimental set-up 

The data sample used in this analysis corresponds to an integrated luminosity C = 114.1± 
2.3 ph-\ collected by the ZEUS detector in the years 1996-2000. During the 1996-97 
data taking, HERA provided collisions between an electroiJl] beam of Ee = 27.5 GeV and 
a proton beam of Ep = 820 GeV, corresponding to a centre-of-mass energy ^/s = 300 GeV 
(£300 = 38.0 ± 0.6 pb~^). In the years 1998-2000, the proton-beam energy was Ep = 
920 GeV, corresponding to = 318 GeV (£318 = 76.1 ± 1.7 pb"^). 

A detailed description of the ZEUS detector can be found elsewhere [13]. A brief outline 
of the components that are most relevant for this analysis is given below. 

Charged particles are tracked in the central tracking detector (CTD) [14], which oper- 
ates in a magnetic field of 1.43 T provided by a thin superconducting coil. The CTD 
consists of 72 cylindrical drift chamber layers, organized in 9 superlayers covering the 
polar-angl^ region 15° < 6* < 164°. The transverse-momentum resolution for full-length 
tracks is (t{pt)/pt = 0.0058pr © 0.0065 © 0.0014/pt, with pr in GeV. 

The high-resolution uranium-scintillator calorimeter (CAL) [15] consists of three parts: 
the forward (FCAL), the barrel (BCAL) and the rear (RCAL) calorimeters. Each part 
is subdivided transversely into towers and longitudinally into one electromagnetic sec- 
tion and either one (in RCAL) or two (in BCAL and FCAL) hadronic sections. The 
smallest subdivision of the calorimeter is called a cell. The CAL energy resolutions, 
as measured under test-beam conditions, are a{E)/E = 0.18/-\/E for electrons and 
a{E)/E = 0.35/ ^/E for hadrons, with E in GeV. 

The position of electrons scattered at small angles to the electron beam direction was 
measured using the small-angle rear tracking detector (SRTD) [16, 17]. The SRTD is 
attached to the front face of the RCAL and consists of two planes of scintillator strips, 
arranged orthogonally. The strips are 1 cm wide and 0.5 cm thick. 

The muon system consists of rear, barrel (R/BMUON) [18] and forward (FMUON) [13] 
tracking detectors. The B/RMUON consists of limited-streamer tube chambers placed 

^ Electrons and positrons are not distinguished in this paper and are both referred to as electrons. 

^ The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
proton beam direction, referred to as the "forward direction" , and the X axis pointing left towards 
the centre of HERA. The coordinate origin is at the nominal interaction point. The pseudorapidity 
is defined as 77 = — In (tan |), where the polar angle, 6, is measured with respect to the proton beam 
direction. 
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behind the BCAL (RCAL), inside and outside a magnetized iron yoke surrounding the 
CAL. These chambers cover polar angles from 34° to 135° and from 135° to 171°, respec- 
tively. 

The luminosity was measured using the bremsstrahlung process ep ejry. The resulting 
small- angle energetic photons were measured by the luminosity monitor [19], a lead- 
scintillator calorimeter placed in the HERA tunnel aX Z = —107 m. 

3 Data Selection 

The data were selected online via a three-level trigger system through a combination of 
three different trigger chains: 

• a muon reaching the inner B/RMUON chambers and matched to a minimum ionizing 
energy deposit in the CAL or any muon reaching the outer B/RMUON chambers or 

• a. D* candidate [20] or 

• a scattered-electron candidate in the CAL [11]. In part of the data taking, the cuts 
on the electron candidate were relaxed if a muon in the inner B/RMUON chambers 
was detected. 

Due to this redundancy, the trigger efficiency for beauty events was high, 94 ± 3% for the 
inclusive study, and 98 ± 2% for the DIS selection. 

Muons were reconstructed offline using the following procedure: a track was found in 
the inner B/RMUON chambers, then a match in position and angle to a CTD track 
was required. In the bottom region of the detector, where there are no inner chambers, 
the outer chambers were used instead. If a match was found to both inner and outer 
chambers, a momentum-matching criterion was added. 

The angular coverage of the B/RMUON and of the track requirements in the CTD restrict 
the muon acceptance to the pseudorapidity region 

- 1.75 <r]^' < 1.3 . (1) 

A cut on the muon transverse momentum 

p^>1.4GeV (2) 

was applied, reflecting the requirement that the muon reaches the inner muon chambers 
in the barrel region. In order to have uniform kinematic acceptance, this cut was also 
applied in the rear region. 
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D* candidates were reconstructed in the D*^ — > > i^r~7r~'")7r+ decay channel (+c.c.) 
making use of the AM (= M{K7nTs) — M[Kii)) technique described in previous pubhca- 
tions [20] with the following cuts: 

mass 1.81 < M{Kn) < 1.92 GeV; 
D* - mass difference 0.1435 < AM < 0.1475 GeV; 
D* transverse momentum > 1.9 GeV; (3) 

D* pseudorapidity \rj^ \ < 1.5; 
K, IT transverse momentum p^'^ > 0.5 GeV; 

slow pion pj= > 0.125 GeV. 

To allow the background to the D* signal to be determined, candidates with wrong- 
charge combinations, in which both tracks forming the candidates have the same 
charge and the third track has the opposite charge, were also retained. 

The hadronic system was reconstructed from the calorimeter information and the re- 
constructed vertex. A four-momentum {p\ , py , p^^ , E"^) was assigned to each calorimeter 
cell. Global hadronic variables were reconstructed by summing over these cells. In the 
case of identified DIS events (see below), the scattered electron candidates were excluded 
from this sum. The inelasticity y was reconstructed from the Jacquet-Blondel estimator 
l/jB = {E — Pz)/2Ee [21], where E — Pz = YIA^^ ~ Pz) tuhs over all cells. In 

the case of DIS events, the alternative value i/e = 1 — ~ cos 6e) as well as the photon 
virtuality were obtained from the energy E'^ and scattering angle 6^ of the final-state 
electron candidate [11]. A sample of events with one muon and one D* candidate was 
selected by requiring: 

• > 1 muon in the muon chamber regions defined by Eq. ([1]) and Eq. ([2]); 

• > 1 D* candidate in the D* acceptance region defined by Eq. Q; 

• the muon candidate track is not one of the three D* candidate tracks, eliminating 
backgrounds from semileptonic D° decays; 

• the D*p system carries a significant fraction of the total transverse energy of the 
event, p!^ ^ / Ej- > 0.14, where Enp is the transverse energy measured by the GAL 
outside a cone of 10° around the proton beam direction to exclude the proton remnant, 
and p^ ^ is the transverse momentum of the D*ii system, reducing combinatorial D* 
background; 

• a reconstructed vertex compatible with the nominal interaction point, suppressing 
non-ep background. 

After this selection, a sample of 232 events remained. The resulting AM distributions for 
the like and unlike D* ^ charge combinations, before the AM cut, are shown in Figs. [2^ 
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and[2]3. 

A subsample of photoproduction events was selected by requiring: 

• no scattered-electron candidate found in the CAL; 

• E-Pz <34: GeV; 

retaining 81% of the inclusive sample. After the unfolding of the detector response, these 
cuts correspond to an effective cut < 1 GeV^ and 0.05 < y < 0.85. The lower limit on 
y arises from the interplay between the 6-quark mass and the acceptance in rapidity. 

Alternatively, a clean DIS sample was obtained by applying the following additional con- 
ditions [11]: 

• a reconstructed electron with energy E'^ >10 GeV; 

• g2 > 2 GeV^; 

• inelasticity yjB > 0.05 and ye < 0.7; 

• AO<E-Pz <60 GeV; 

• the electron hits the rear calorimeter outside a rectangle of |X| < 13 cm and |y| < 7 
cm. 

These cuts correspond to an effective cut > 2 GeV^ and 0.05 < y < 0.7. For this 
sample, which contains less combinatorial background, the D* cuts were relaxed to 

1.5 GeV; 

0.4 GeV; (4) 
0.12 GeV; 

and the cut on p!j^ ^ / Et was dropped. All other cuts on the D* and the muon remained 
unchanged. A sample of 44 events was obtained. The resulting AM distributions for the 
like and unlike D* n charge combinations are shown in Figs. Eb and[2li. 

4 Backgrounds and event simulation 

Several contributions to the selected data sample were evaluated: 

• the signal from beauty decays; 

• the background from fake D* combinations; 

• the D* ^ background from charm decays; 

• the background from fake or non-prompt muons with a real D* from charm. 



D* 

Pt 
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For the signal from beauty and charm production, Monte Carlo (MC) simulations were 
performed using the Pythia [22], Rapgap [23] and Herwig [24] generators. These 
simulations are based on leading-order matrix elements complemented by parton showers 
to obtain higher-order topologies. The direct photon-gluon fusion process (7(7 — QQ., 
Q = b,c), flavour excitation in the resolved photon and proton (e.g. Qg Qg, '-/Q Qg), 
and hadron-like resolved photon processes {gg QQ) were included. Gluon splitting into 
heavy flavours {g — > QQ) in events with only light quarks in the hard scattering was not 
included in the simulations; this contribution is, however, expected to be small [25]. For all 
generated events, the ZEUS detector response was simulated in detail using a programme 
based on GEANT 3.21 [26]. 

The number of background events under the D* mass peaks (fake D*) was estimated using 
the wrong-charge Ktt combinations combining the like- and unlike-sign samples. This was 
found to minimize the bias due to charge correlations [25]. Wrong-charge combinations 
were normalised to the data outside the D* peak in the side-band 0.15 < AM < 0.17 GeV, 
separately for the like-sign and unlike-sign D*fi sample, as shown in Fig. [2j Dedicated 
studies [25] performed by selecting data on the D* side-band showed that this procedure 
correctly reproduces shape and normalisation of the iak.e-D* background for the relevant 
variables used in the analysis. 

Fake muons can be produced by hadron showers leaking from the back of the calorimeter 
or by charged hadrons traversing the entire calorimeter without interaction. In addition, 
low-momentum muons can originate from in-fiight decays of pions and kaons. It is also 
possible for tracks reconstructed in the central tracker to be wrongly associated to a signal 
from a real muon in the muon chambers. A dedicated study [25] based on pions from K'^ 
decays, protons from A decays, and kaons from and D* decays, showed that the detector 
simulation reproduced these backgrounds reasonably well. The fake muon probability for 
the — > TT+vr" sample is about 0.2%. Most fake muons are associated with fake D* 
candidates, and therefore accounted for in the {ak.e-D* background estimated directly 
from the data. Fake muons associated with a real D* are included in the charm and 
beauty MC samples. 

Distributions of variables used in the event selection or relevant for the event kinematics, 
for the unlike-sign inclusive sample, are compared to the expectations from these simula- 
tions in Fig. [3], separately for the beauty- and charm-enriched regions defined in Section 
[5l Agreement with expectations is obtained, apart from some possible deviations in the 
^ and '^/Et distributions in the beauty-enriched region, which are accounted for in 
the systematic uncertainties (Section [7]). 
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5 Signal extraction 



In this section, the signal-extraction procedure is described for the inclusive sample. The 
7P subsample and the DIS sample were treated in an analogous way. 

Figures and show the distribution of the angular difference AR = -\/A0^ + Ar/^ 
between the D* and the muon, where (p is the azimuthal angle, for events passing all 
selections, including the AM cut. The distributions are shown separately for like- and 
unlike-sign D*fi events. The expected signal and background distributions, normalised 
to the fractions determined later in the analysis and described below, are also indicated. 
For unlike-sign events, the region AR > 2, which mainly corresponds to the back-to- 
back configuration, is clearly dominated by events from charm. Indications that the 
simulated distribution might be more sharply peaked than the data turned out to have 
little influence on the determination of the beauty fraction. In contrast, the region AR < 2 
is enriched in beauty events, in which the D* and muon originate mainly from decays of the 
same parent B hadron. This is illustrated further in the D*fi invariant-mass distribution 
(Figs. and HJi) for events in the beauty-enriched region {AR < 2). A peak with an 
upper edge close to 5 GeV, which can be attributed to the partial reconstruction of the 
decaying B meson, is clearly visible. A comparison with the like-sign sample shows that 
the low-mass edge of this peak is dominated by background. An invariant-mass cut of 
3 GeV < M{D*fi) < 5 GeV was therefore applied to the AR < 2 subsample. 

After this additional cut, and after statistical subtraction of the iake-D* background, the 
contributions of charm and beauty were determined by a two-component fit to the AR 
or A(f) distributions, shown in Figs. [5^ and[5l3. The fake- muon background with a real D* 
from charm, which is treated as part of the charm component, is shown separately. The 
small fraction of fake muons from beauty was included in the beauty component. The fit 
result for the fraction of beauty events in the final inclusive sample shown in Figs. [5^ and 
[Sb, using the shapes predicted by Pythia, is: 

• fb = 0.307 ± 0.064(stat.) for the AR and 

• fb = 0.290 ± 0.062(stat.) for the A0 distribution. 

The AR result was chosen as the reference, and the A0 result used as a systematic 
check. With these fitted fractions, the breakdown into the corresponding number of 
beauty, charm, and fake-muon candidates in each subsample is given in Table [TJ In the 
unlike-sign part, the beauty and charm contributions are well separated, with only small 
cross-contaminations. The normalisation of the beauty and charm contributions in Fig. H] 
already reflects these fitted fractions. Agreement is seen, also in the like-sign part, which 
was not included in the fit. 

The results from the application of the same procedure to the 7p subsample are also 
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shown in Table [H The analogous results for the DIS sample are shown in Figs. \5jp and 
|5ll. The A(f) distribution gives less discrimination in this case, due to the transverse boost 
from the exchanged virtual photon. Therefore, the M{D*fi) distribution was used. The 
fitted beauty fractions in the DIS sample, using the shapes predicted by Rapgap, are 

• fb = 0.55 ± 0.25(stat.) for the AR and 

• fb = 0.43 ± 0.30(stat.) for the M{D*fi) distribution. 

Again, the AR result is chosen as the reference, and the other as a cross check. The 
breakdown into different event categories is shown in Table [TJ The acceptance corrections 
for the cross sections which will be presented in Section [8] were evaluated using Pythia 
for < 1 GeV^, Rapgap for > 1 GeV^, and Herwig as a systematic check. 

6 Theoretical predictions and uncertainties 

For direct comparisons with QCD, next-to-leading-order predictions were used. Calcula- 
tions in which b quarks are treated as massless particles [27] are not applicable in this 
kinematic range, while calculations based on alternative parton-evolution schemes [28] do 
not yet exist with full NLO implementation. Fixed-order NLO calculations with mas- 
sive b quarks should yield accurate predictions. Different types of such predictions were 
evaluated. 

The FMNR program [1] evaluates cross sections for next-to-leading-order beauty pro- 
duction in 7p collisions in the fixed-order massive approach, for both point-like and 
hadron-like photon coupling to the heavy quarks. The parton-density functions used were 
CTEQ5M [29] for the proton and GRV-G-HO [30] for the photon. The renormalisation 
and factorisation scales n were chosen to be equal and parametrised by /io = \/pt + 
where p\ is the average of the squared transverse momentum of the two emerging b quarks 
and TTib = 4.75 GeV. An estimate of the theoretical uncertainty was obtained by simulta- 
neously varying 4.5 < nii, < 5.0 GeV and /io/2 < /i < 2/io such that the uncertainty was 
maximised. Typical uncertainties resulting from this procedure (e.g. for the bb total cross 
section) are +40% and -25%. Variations of the parton densities led to uncertainties which 
were much smaller than the uncertainties related to mass and scale variations. They were 
therefore neglected. 

Predictions at the level of visible D*fi final states are needed in addition to those at 
parton level. The FMNR program provides a framework to fragment b quarks into B 
hadrons, and simulate the decay of these hadrons by interfacing them to appropriately 
chosen decay spectra. However, decays to complex final states, such as a D* and fi from 
the same B hadron with cuts on both particles, cannot be easily implemented in this 
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scheme. A straightforward interface of the parton-level events produced by FMNR to 
MC-hke fragmentation and decay chains is also impracticable, since these events have 
positive and negative weights spanning more than 8 orders of magnitude, making such an 
approach extremely inefficient. 

These difficulties were overcome in a two-step process. In the first step, two or more 
FMNR parton-level events with large positive and negative weights and similar topology 
were combined into events with much smaller weights by averaging the parton momen- 
tum vectors [31]. Events were considered to have similar topology if the differences in 
transverse momentum, rapidity and azimuthal angle of the b quarks were less than user 
cut values that reflect the detector resolution. Furthermore, events with small weights 
were sampled with a probability proportional to their weight. In this way, the weight 
range was reduced to about two orders of magnitude. It was explicitly checked that this 
procedure preserves the NLO accuracy for the relevant cross sections at parton level (e.g. 
b quark pt and angular distributions). 

In the second step, these parton-level events were interfaced to the Pythia/Jetset [32] 
fragmentation and decay chain, making use of the full decay tables and decay kinematics 
implemented in Pythia 6.2. Therefore, non-dominant complex decays, such as B ^ 
D*D followed by ^ /xX, or muons produced through intermediate J/ip oi t states, 
were automatically included. The initial-state partons were allowed to have intrinsic kx 
(typically ~ 300 MeV) as implemented in Pythia. This has a negligible effect on the 
resulting cross sections (~ 1%). Parton showering was not included in order to avoid 
double counting of higher-order contribution^. 

Fragmentation of b quarks close to production threshold is non-trivial. The details of 
the threshold treatment were found to be much more important than the choice of a 
particular fragmentation function. The Peterson formula [34] with e = 0.0035 was used. 
Three approaches were considered: 

• independent fragmentation as implemented in Pythia [32]. The use of this quite old 
model was motivated by the fact that FMNR does not provide colour connections on 
an event-by-event basis; 

• fragmentation in the Lund string model [35], again as implemented in Pythia. For 
this purpose, reasonable colour connections were assigned to each event; 

• the independent fragmentation scheme provided within the FMNR framework, rescal- 
ing the 5-hadron momentum to a fraction of the 6-quark momentum according to the 
Peterson formula, which is a somewhat crude approximation at threshold. 

^ The MC@NLO approach [33], which aUows the combination of NLO matrix elements with parton 
showers, is not yet available for ep interactions. 
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The second option was used for all central predictions. The first option was used to obtain 
the lower systematic error (typically -5%). The third option could not be applied directly, 
since it does not provide cross-section predictions for correlated final states from the same 
b quark, as needed here. Instead, it was applied to a cross section in which the final-state 
correlations originate from different b quarks only, which is more easily calculable in this 
scheme. The results were used to evaluate a generic upper systematic error of +15% on the 
fragmentation procedure close to b production threshold. The effect of a variation of the 
Peterson parameter e in the range 0.0023 to 0.0045 was found to yield uncertainties that 
were much smaller than the uncertainties due to the different fragmentation procedures. 
It was therefore neglected. 

The branching fractions were corrected to correspond to those obtained from the Parti- 
cle Data Group [36], as listed in Table [2l Branching fraction uncertainties resulted in 
uncertainties on the D*fi cross section of typically ±12%. 

In principle, FMNR predictions are only valid for the photoproduction regime. The 
Weizsacker- Williams approximation with an effective Q^ax cutoff of 25 GeV^ (~ ml) [37] 
was used to include the ~15% DIS contribution to the combined cross section. 

Alternatively, the DIS part can be calculated using the NLO predictions from HVQDIS [2]. 
Only point-like contributions are included in this prediction. The parton density function 
used was CTEQ5F4 [29] . The renormalisation and factorisation scales n were chosen to be 
equal and parametrised by /xq = \/Q'^ + '^1- Mass and scales were varied as for FMNR. A 
scheme for the calculation of visible cross sections for correlated final states, corresponding 
to the FMNR®Pythia interface described above, was not available. Therefore, cross- 
section comparisons in DIS are limited to the parton level. 

7 Systematic uncertainties 

The main experimental uncertainties are described below, in order of importance. Num- 
bers in parentheses are quoted for the inclusive selection. Uncertainties for the 7p results 
are also quoted when they differ significantly from the inclusive results. For the DIS 
sample, the statistics were often too small to derive meaningful systematic errors. The 
errors from the inclusive sample were used instead. 

• Simulation of ^ . The largest error arises from the observation that the muon and 
D* pt spectra in the b signal region of the data {AR{D*n) < 2, 3 < M{D*n) < 5 GeV) 
appear to be somewhat softer than predicted by the Monte Carlo simulations (Fig. [3]). 
The differences are concentrated at small values of pf^ ^ / Et- Since the corresponding 
spectra are well reproduced in the charm region with larger statistics, this cannot be 
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attributed to problems with the muon or D* reconstruction. There are several ways 
to interpret these differences: 

a) they are statistical ffuctuations. This assumption leads to the central result re- 
ported; 

b) the signal distribution is significantly softer than predicted by QCD. Due to the 
rising efficiency as a function of p^, this would change the efficiency calculation for 
the measurement of the visible cross section. To evaluate this possibility, the MC 
p!^ ^ (true level) distribution in the signal region was reweighted to be compatible 
at the 1 a level with the measured ^ spectrum of the inclusive sample (Fig. ^ 
(+14%); 

c) there is an additional unknown background contribution at low pj-, which occurs 
only in the beauty-enriched region. There is no indication that this is the case. 
Nevertheless, to account for this possibility, the ^/ Et cut was tightened from 
0.14 to 0.2, which removes most of the differences (Fig. [3]) (-33% for inclusive, 
-18% for 7p selection). 

Branching fractions. The beauty-enriched region, which dominates the fit result, 
is mainly populated by events in which the D* and fi originate from the same b. The 
rate of these events depends on different branching fractions from those relevant to 
the charm-enriched region, in which the D* and fi originate from different b quarks. 
A variation of these branching fractions, within the uncertainties quoted in Table Ej 
therefore affects the shape of the beauty contribution and the fitted beauty fraction 
(±8%). 

Fragmentation and parton showering. The Herwig MC uses a different fragmen- 
tation model from that of Pythia and Rapgap. It also yields different bb correlations 
from direct /resolved contributions and parton showering. This leads to differences in 
the acceptance, and in the fitted beauty fraction (+5/-8%). 

Signal-extraction procedure. In addition to statistical fluctuations, different ways 
to fit the data can yield systematic differences due to binning effects and different 
systematics for different variables, e.g. imperfections in the shape of the MC distri- 
butions. To check the error from this effect, the cross sections were evaluated using 
different procedures: fits to AR{D*fi), A0(D*/i), M{D*fi), and simple event count- 
ing. In all cases the differences were well within the quoted errors. To avoid double 
counting of statistical and systematic errors, these were used as cross checks only. 

Uncertainty on the estimation of the muon chamber efficiency. Corrections 
to the MC muon chamber reconstruction efficiency were obtained from independent 
data samples and varied within their uncertainties (±5%). 
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• Fake muon background. The background from fake muons has been extensively 
studied [25] and is further constrained by the hke-sign distribution of Fig. |Dd, which 
is dominated by this background. Accordingly, it was varied by a factor 1.5 (-4%). 

• Luminosity measurement. The uncertainty associated with the luminosity mea- 
surement for the 1996-00 data taking periods used for this analysis was included 
(±2%). 

• Tracking. All tracking-based cuts (pr and mass cuts) were varied by their respective 
uncertainties. To avoid double counting of statistical uncertainties, the D*-related 
systematics were taken from previous ZEUS DIS [38] and 7p [39] analyses employing 
similar cuts but with larger event samples. The cut on was varied by ib40 MeV. 
This yielded a combined error of +6% and -4%. 

• Trigger acceptance. The error on the trigger acceptance was evaluated by comparing 
the efficiencies of the different trigger chains in the data with each other and with the 
MC (±3%). 

• B^-B^ mixing. The possible systematic effect due to the variation of the mixing rate 
was found to be negligible. 

The total systematic uncertainty was obtained by adding the above contributions in 
quadrature. 

8 Results 

To present results from the combined data sets, the measurements from the 1996-97 run 
at = 300 GeV have been corrected using the predicted cross section ratio [1] of 1.06, to 
correspond to the higher center-of-mass energy of 318 GeV. All cross sections are therefore 
quoted for ^/s = 318 GeV. 

8.1 Visible cross sections 

The first step is the extraction of visible cross sections for the D*fi final state from beauty. 
The acceptance for the D* D^tTs — ^ {K'iT)'n's decay chain was unfolded using a branching 
fraction of 2.57 ±0.06% [36]. The effective b branching fractions used in the different MC 
generators were corrected to those listed in Table [2] in order to account for their influence 
on the overall acceptance, and on the shape of the predicted beauty contributions. 

The measured beauty fraction in the inclusive sample, corrected for detector acceptance 
and branching fractions, was used to obtain the cross section for the process ep — ebbX — > 
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eD*^^X in the visible kinematic range p^* > 1.9 GeV, —1.5 < 7]^ < 1.5, > 1.4 GeV 
and —1.75 < r]'^ < 1.3 as: 

(Tvisiep ebbX eD*^fxX) = 160 ± 37(stat.)+5°(syst.) pb. (5) 

This includes both unlike- and like-sign D*fi combinations. The leading-order cross sec- 
tions predicted by Pythia and Herwig in the same kinematic range are o"vis(ep — >■ 
ebbX — >■ eD*nX) = 80 and 38 pb, respectively. The measured cross section is larger than, 
but compatible with, the FMNR®Pythia NLO prediction 

^vi^°(ep ^ ebbX ^ eD*^fxX) = 67+?;(NLO)+;'(frag. © br.) pb, (6) 

where the first error refers to the uncertainties of the FMNR parton-level calculation, and 
the second error refers to the uncertainties related to fragmentation and decay. 

For the photoproduction subsample, a visible cross section in the kinematic range < 1 
GeV^ and 0.05 < y < 0.85 was obtained: 

cTvis,7p(ep ebbX eD*^fiX) = 115 ± 29(stat.)+2j(syst.) pb. (7) 

This can be compared to the NLO prediction from FMNR®Pythia, 

^S?p(eP e^^^ ^ eD*^fiX) = 54l}g(NLO)+;°(frag. © br.) pb. (8) 

As in the inclusive case, the NLO prediction underestimates the measured cross section 
by about a factor of 2, but is compatible with the measurement (Table [3]). 

From the DIS sample, a visible cross section in the kinematic range > 2 GeV^, 0.05 < 
y < 0.7 and p^* > 1.5 GeV (other D* and muon cuts as for Eq. ([5])) of 

c^vis,Dis(ep ^ ebbX eD*fxX) = 58 ± 29(stat.)+2j(syst.) pb (9) 

was obtained. 

Again, the cross sections obtained from Rapgap (used to compute acceptance corrections 
for the central signal extraction) and Herwig (used for systematic checks, particularly 
with regard to differences in the bb correlations) in the same kinematic regime are consid- 
erably lower, a{ep ebbX eD*fiX) = 26 and 10 pb, respectively. An NLO prediction 
is not available for this kinematic region. 

8.2 Comparison to HI results 

A photoproduction cross section similar to Eq. ([7]) in a slightly different kinematic range, 
p!j^* > 1.5 GeV, \ri^'\ < 1.5, p'^ > 2.0 GeV, \r]''\ < 1.735, < i GeV^ and 0.05 <y< 0.75 
has been obtained by the HI collaboration [12] 

c^vis,7p(ep ebbX eD*^fiX; HI) = 206 ± 53(stat.) ± 35(syst.) pb. (10) 
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The ZEUS cross section of Eq. ([7]) extrapolated to the same kinematic range as the HI 
measurement using FMNR(8)Pythia is 

(Xvis,7p(ep ebbX ^ eD*^fiX; HI) = 135 ± 33(stat.)^3t(syst.) pb, (11) 

which is somewhat smaller, but in agreement within errors. 
The corresponding FMNR®Pythia NLO prediction is 

t^vis,?p(ep ebbX ^ eD*^fiX; HI) = 6ltJ^(NL0)+f (frag. © br.) pb. (12) 

This is larger than the NLO cross section evaluated by HI [12] due to the inclusion of the 
hadron-like photon contribution, the inclusion of secondary-muon branching fractions for 
D* and ^ from the same b quark (Table [2]), and a detailed simulation of the kinematics 
of the b ^ B ^ D* chain rather than direct coUinear fragmentation of b quarks into 
D* mesons. The data to NLO ratio is consistent with the results in the ZEUS kinematic 
range. 

8.3 Cross sections for D*fi from the same b quark 

In all the cross sections evaluated above, a significant part of the systematic error arises 
from the fraction of the beauty contribution in the charm enriched or like-sign regions, 
where it cannot be well measured {AR > 2 region in Fig. H^; Figs. Hb andHJi). This 
fraction depends on details of the description of bb correlations in the MC used for the 
signal extraction. In the beauty-enriched low-Ai? unlike-sign region, which dominates the 
fit of the beauty fraction, about 95% of the D*fi pairs are produced from the same parent 
b quark. The systematic error can thus be reduced by reinterpreting the measurements 
in terms of cross sections for this subprocess only. The corresponding cross section for 
photoproduction of a D* and n from the same b quark (always unlike sign, same kinematic 
cuts as for Eq. ([7])) is 

crvis,7p(ep ^ ebib)X, b{b) eD*fiX) = 52 ± 13(stat.)+Ji(syst.) pb (13) 

where b{b) stands for the sum of b and b cross sections, and all other cuts remain the 
same. 

This can be compared with the NLO prediction 

^ eb{b)X, b{b) ^ eD*fiX) = 29l^(NLO)+;(frag. © br.) pb. (14) 

For the DIS kinematic range (same as Eq. iQ) 

c^vis,Dis(ep ^ eb{b)X, bib) eD*fxX) = 28 ± 14(stat.)+^o(syst.) pb (15) 
is obtained. 
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8.4 Parton- level cross sections 



For a direct comparison with the NLO parton-level predictions, the measured visible cross 
sections were extrapolated to 6-quark level. In order to minimize the systematic error, the 
6-level cross section is quoted for individual b (or b) production rather than for correlated 
bb-paiT production, i.e. using the cross sections displayed in Eqs. (1131) and (fT5|) . 

A significant fraction of the parent b quarks of the selected events is expected to have very 
low values [25]. Therefore, cross sections with no cut on have been measured. Fur- 
thermore, there is a strong correlation between the pseudorapidity of the D*n system and 
the rapidity of the parent b quark, = ^^3^. In order to reflect the limited angular 
acceptance of the detector for both the D* and the muon, the cross-section measurement 
was restricted to < 1. In this range, restricted to photoproduction, the p\, and distri- 
butions of Pythia (after parton showering) agree with the central NLO 6-quark spectra 
from FMNR to within ±15% [25]. Therefore, Pythia was used to extrapolate the visi- 
ble cross section for the photoproduction region. Similarly, the corresponding Rapgap 
spectra for the DIS case agree [25] with the central NLO predictions from HVQDIS. 

The acceptance for b quarks due to the kinematic cuts on the fragmentation and decay 
products ranges from ~ 4% at = GeV to ~ 55% at p^ = 10 GeV. The remaining 
part of the extrapolation is due to the relevant branching ratios. 

The extrapolation implies additional systematic uncertainties from the fe-quark fragmen- 
tation (+5/-15%) and decay (±9%) and the details of the shape of the p^ spectrum 
(±5%). The extrapolation was calculated assuming the validity of the NLO p^ shape and 
is therefore valid only in the context of this theoretical framework; the uncertainty for 
the visible cross section corresponding to a potential deviation from this shape, namely 
the reweighting of the ^ spectrum, is removed. The result for the extrapolated cross 
section for < 1, < 1 GeV^ 0.05 <y< 0.85 and = 4.75 GeV was 

a^p{ep b{b)X) = 11.9 ± 2.9(stat.)+J;3(syst.) nb. (16) 

The corresponding result for the extrapolated cross section for < 1, > 2 GeV^, 0.05 < 
y < 0.7 and mf, = 4.75 GeV was 

(^Disiep b{b)X) = 3.6 ± 1.8(stat.)+J;5(syst.) nb. (17) 

These cross sections are to be compared to the NLO prediction for the same kinematic 
range using the FMNR calculation of 

</«(ep ^ b{b)X) = 5.8tl-l nb, (18) 
and to the NLO HVQDIS prediction of 

^DiFiep - b{b)X) = 0.87tUl nb. (19) 
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These cross sections are presented in Fig. O The ratio of measured to predicted cross 
sections in the photoproduction region remains the same as the ones obtained from the 
comparison at visible level (Table [S]). This confirms the self-consistency of the extrapola- 
tion procedure used. 



8.5 Comparison to previous ZEUS measurements 

In order to compare the photoproduction cross section to previous ZEUS results [8,10], 
the cross sections Eq. f[T5]) or equivalently Eq. flTB]) . already referring to the production 
of a single b quark, need to be translated into a differential cross section, in the 
pseudorapidity range \ri^\ < 2 [10]. The median p\, value for events satisfying the cuts for 
Eq. (fT5l) is 6.5 GeV [25]. The measured cross section, Eq. (ITHIl . is therefore extrapolated 
to this value using FMNR(g)PYTHlA, yielding 

-nr{PT = 6-5 GeV, \Tf\ < 2) = 0.30 ± 0.07(stat.)+J;Jg(syst.) nb. (20) 

(XDrp 

This result is compared to theory and previous measurements in Fig. [71 It is higher than, 
but consistent with, these measurements. 



9 Conclusions 



Cross sections for beauty production in ep collisions at HERA have been measured in both 
the photoproduction and DIS regimes using an analysis technique based on the detection 
of a muon and D* . Agreement is obtained with the corresponding HI result. Since the 
analysis is sensitive to 6-quark production near the kinematic threshold, the measured 
visible cross sections were extrapolated to 6-quark cross sections without an explicit cut 
on p\,. Both at visible and at quark level, the measured cross sections exceed the NLO 
QCD predictions, but are compatible within the errors. The data to NLO ratio is also 
larger than, but compatible with, previous ZEUS measurements of the h production cross 
section at higher p\. 
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Sample 








Charm 




Cuts 


Data 


Beauty 




fake /i 


Fake D* 






prompt yU 






Inclusive sample 


unlike sign 


Ai? < 2, 3 < M < 5 GeV 


41 


28.8 


3.6 


2.5 


11.2 
















Ai? > 2 


93 


6.5 


56.5 


18.4 


6.4 


like sign 


Ai? < 2, 3 < Af < 5 GeV 


18 


1.6 


0.7 


1.3 


9.2 
















Ai? > 2 


36 


11.1 


2.0 


21.2 


5.2 




IP 


unlike sign 


Ai? < 2, 3 < M < 5 GeV 


31 


23.6 


1.2 


1.6 


8.3 
















Ai? > 2 


79 


6.2 


48.8 


14.0 


6.2 


like sign 


Ai? < 2, 3 < M < 5 GeV 


14 


1.5 


0.6 


0.4 


6.9 
















Ai? > 2 


28 


9.1 


1.9 


17.1 


5.2 




DIS 


unlike sign 


Ai? < 2, 3 < M < 5 GeV 


11 


8.1 


1.6 


0.5 


1.0 
















Ai? > 2 


14 


3.2 


6.1 


1.2 


3.9 


like sign 


Ai? < 2, 3 < M < 5 GeV 


3 


1.2 


0.1 


0.5 


0.3 
















Ai? > 2 


10 


5.0 


0.2 


2.6 


1.2 



Table 1: Composition of the final D*fi event samples (number of events) as 
determined from the fit to the AR distribution. 
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O 1^ in /~\ 1 

CllcLIlIlcl 


l^T* O TO /^To 1 TO *T "fv* O T TO 7" / 7-^ ^ 7-^ TOO T X/'T TO 

Diancnmg iiaciion w/o lj -d niixuig 


6 ^ D*+ inclusive 


i ( .O It z.u /o 

OD ± O /o iJ , 14 ± O /o iJ 


6 — j> jj, direct 


iu.yo It u.z ( /o 


Ji V //+ 1 n r~l 1 "vo r*"!" 

U ^ Li lliQlicCL 


97 _|_ f) 4n 0/ 

o.z 1 HI u.^u /o 


b fi~ indirect 


2.21 ± 0.50 % 


all b fi^ 


21.43 ±0.70 % 


bb D*^iJ,^ (diff. 6's) 


4.34 ± 0.92 % 


bb D*^iJ,^ (diff. 6's) 


3.08 ± 0.60 % 


b direct 


2.75 ±0.19 % 


6 D*^^"^ indirect 


1.09 ±0.27 % 


all b D*^^^ 


3.84 ± 0.33 % 



Table 2: Branching fractions assumed for cross-section determinations. The 
indirect contributions include cascade decays into muons via charm, anticharm, r+ 
and J/ip. The values in the Table are given before the inclusion of the effect of 
B^-B^ mixing (mixing parameter x = 0.1257 ± 0.0042 ) [36]. 



cross section 


measured 


NLO QCD 


ratio 


inclusive, visible, Eq.(l5]) 
7P, visible, Eq.(j7]) 
DIS, visible, Eq.© 


160 ± 37tl7 pb 
115 ±29^^^ pb 
58 ± 29tll pb 


Q7+_it pb 
54tl^ pb 


2 4+0.9 

1 +0.8 
^•-■--i.o 


7P, vis. same b, Eq.(fT3l) 
DIS, vis. same b, Eq.(fT5!) 


52 ± 13tu pb 
28 ± 141^0 pb 


29tl pb 


1 0+0.7 

^■^-0.8 


7P, b quark, Eq.(fT6|) 
DIS, b quark, Eq.([T7D 


11.9±2.9l^j nb 
3.6± 1.8l?i nb 


S-S^ii nb 
0.871°:?^ nb 


2.0+?:? 


7P, b differential, Eq.(|20|) 


0.30 ± 0.07tom nb/GeV 


0.16^°;°^ nb/GeV 


■'■•«-0.8 



Table 3: Comparison of measured and predicted cross sections. For the measured 
cross sections, the first error is statistical, and the second systematic. For the QCD 
prediction, the error is due to the parton-level NLO calculation convoluted with the 
uncertainties of fragmentation and decay to the visible final state. The number in 
parentheses refers to the corresponding equation for each cross section, see text. 
For the definition of the kinematic range of each cross section see text. 
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Figure 1: Processes leading to D* ji final states. 
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ZEUS 




0.14 0.15 0.16 0.17 0.14 0.15 0.16 0.17 



M(K7t7t3)-M(K7t) (GeV) M(K7t7t3)-M(K7t) (GeV) 

Figure 2: Distribution of AM for data (full circles) and combinatorial background 
(hatched histogram) for (a) inclusive unlike-sign (b) inclusive like-sign (c) DIS 
unlike-sign and (d) DIS like-sign muon-D* combinations. The D* signal region is 
indicated by the shaded area. 
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Beauty enriched region 



ZEUS 



Cliarm enriclied region 




Figure 3: Distributions of variables used in the event selection or relevant for the 
event kinematics, for the unlike-sign inclusive sample after all cuts. All variables 
are defined in the text except ntrks (total number of tracks) and C^*^^ (rapidity 
of the D*n system). The distributions are shown separately for the beauty- and 
charm- enriched regions as defined in Section\^ The beauty, charm, fake-muon and 
wrong- charge Kir background distributions are indicated by different shading styles, 
and normalized to the fractions determined later in the analysis. 
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• ZEUS (114 pb"') 

□ charm (prompt n) 
■ beauty 

□ charm (fake [j.) 
^ fake D* 



unlike-sign D {x 



(a)i > 




3 4 

AR(DV) 



5 6 7 

M(dV) (GeV) 



Figure 4: Distribution of (a),(h) AR{D*fi) and (c),(d) M{D*fi) for inclusive 
data (full circles), beauty and charm signal, fake-fi and fake- D* backgrounds. The 
latter are distinguished by different shading styles. Unlike-sign and like-sign D*fi 
combinations are shown separately. Cuts described in the text are indicated by the 
arrows. The relative contributions of charm and beauty are determined by the fit 
to the AR distribution. 



25 



ZEUS 











re 
;g 


12 


■B 




c 
re 


10 


u 






8 




6 




4 




2 









DIS 



3 4 

AR(DV) 



(b) 12 
10 
8 
6 
4 
2 




• ZEUS (114 pb"") (c); 
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Figure 5: Distributions of (a) AR{D*fi), (c) A(f){D*fi) for unlike-sign events in 
the inclusive sample, (h) AR{D*fi) and (d) M{D*n) for unlike-sign events in the 
DIS sample after subtraction of the fake-D* background. Data points (full circles) 
are shown together with the contributions from beauty and charm, as determined 
in the fit. The fake-muon contribution from charm is shown separately, but fitted 
together with the charm contribution. 
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Figure 6: Cross section for single b or b- quark production in the rapidity range < 
1 for photoproduction (left) and DIS (right), compared to NLO QCD predictions 
from FMNR (left) and HVQDIS (right). The ^p cross section is for 0.05 < y < 
0.85, g2 < I QeY2 ^/^g j^jg ^^^gg section for 0.05 <y < 0.7, > 2 GeV'^ . 
The cross sections are integrated over the ranges, and over the full p\ range. 
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Figure 7: Differential cross-section ^ of this analysis (filled square) compared to 

measurements from previous ZEUS analyses, after rescaling to the kinematic range 
indicated in the Figure. 
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